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1. Introduction 

In Lattice QCD, the Wilson discretization of the fermion action introduces 0(a) cutoff effects, 
which may be systematically reduced by the well-established Symanzik improvement programme 
[1]. In the case of massless fermions, only the clover term with the Sheikholeslami-Wohlert coeffi¬ 
cient c sw [2] needs to be added to the action, but composite fields require additional improvement 
terms [3]. 

The function c sv/ (g^) has been determined for various numbers of massless, dynamical quark 
flavours and different actions [4, 5, 6, 7, 8, 9, 10]. It has found successful application in simulations 
with massless quarks (e.g. [11] with Nf = 4) or light quarks (e.g. [12, 13] with Nf = 3). In the 
latter case, in addition to the clover term there exist mass-dependent improvement terms of order 
0(am q light), which can usually be treated perturbatively or even be neglected if am q u g ht is small. 

The present work aims at improvement of the general, non-degenerate four flavor Wilson 
fermion action with a massive charm quark. In that case, an incomplete knowledge of the coef¬ 
ficients of the mass-dependent improvement terms introduces large cutoff effects, see below. We 
want to avoid both these large effects and the non-perturbative determination of a large number of 
improvement coefficients. This is possible by switching to a massive improvement and renormal¬ 
ization scheme, in which the improvement condition for c sw is formulated at finite charm quark 
mass along a line of constant physics (LCP). 

2. Symanzik improvement with four non-degenerate quark flavors 

We are interested in lattice QCD simulations with four flavors of Wilson quarks i//, 

Jz? = Jzf(7 + = w(Dw + 4/ + m cr jt)v^ ) (2.1) 

=^/,o = --2? + aff sw , =^sw = Cswfeo) 5*^4 vFp v yt , (2.2) 

and consider a non-degenerate matrix of bare subtracted quark masses 

M = diag (m qtU ,m q 4 ,m q ji,m qtC ) , m qj = m,, 0 - m crit . (2.3) 

The particular choice of a lattice gauge action Jz?g is irrelevant for the present discussion, but will 
be the tree-level improved Liischer-Weisz gauge action in section 5. While 2zf/ o is the correct 
on-shell 0(a) improved action density for simulations at the critical line (m qq = 0), additional 
mass-dependent 0(a) terms arise for non-vanishing quark masses. They read 

= (7! xjfM 2 xi/ + G 2 Tr(M)\jfM\i/+ (a 3 Tr(M 2 ) + a 4 Tr(M) 2 ) ynjr + a 5 Tv(M)Tv(F^ v ) (2.4) 

and should be included in the simulation (S^i.m = -£/,o + with suitably chosen coefficients 

cjf gjf. One may absorb the above improvement terms into a redefinition of the bare parameters of 
the theory, 

sl go = gl( 1 + ab g (gl)Tx[M]/Nf) , m qJ -> m q f = m qJ fi(gl, {aM}) , (2.5) 

where / is a non-trivial function of the bare subtracted quark masses in aM [3, 14] and a number 
of improvement coefficients related to those in (2.4). 
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For full 0(a) improvement of the action at finite quark mass, a non-perturbative determination 
of the numerous improvement coefficients is in principle necessary. However, this is clearly out 
of reach for any number of flavors. In the past, without a charm quark, one could circumvent 
the problem due to the smallness of the mass-dependent improvement terms. For typical lattice 
spacings a < 0.1 fm, assuming m q s ~ 100 MeV and m qc ~ 1 GeV, they are roughly of sizes 

am q, s < 0.05 , am q c < 0.5 . (2.6) 


Hence, a viable way for small lattice spacings and quark masses, i.e. simulations with Nf < 3, has 
been to treat the terms appearing in as small and account for them approximately. 

Regarding (2.6), the effects of the charm mass in the improved theory are numerically much 
more important, which makes the non-perturbative treatment of all improvement coefficients un¬ 
avoidable. Even if such a treatment was feasible, any uncertainties in the improvement coefficients 
would be amplified by the large charm mass, leading to a significant uncertainty in the improved 
theory. Note that the same issue concerns the renormalization constants, for which it arises already 
at the O(a 0 ) level, 


WlR,i — Z m 


m q ,i + (r m 


Tr[M]' 
j Nf . 


+ O(a) . 


(2.7) 


In particular for observables in the light sector (no charm valence quarks involved), the large charm 
mass cutoff and renormalization effects and assigned uncertainties would be likely to surpass the 
small physical effects of a massive dynamical charm quark [15] that one is interested in. 


3. Massive renormalization and improvement scheme 

For the reasons discussed in the previous section, the massless renormalization and improve¬ 
ment scheme is impractical and unstable if a physical charm quark is to be included in the sea. 
Instead, we choose a massive scheme, in which the mass dependent improvement terms are ab¬ 
sorbed in the bare parameters Eq. (2.5), which are tuned such that a LCP is followed. Accordingly, 
the renormalization and improvement conditions are imposed at finite quark masses, and the coef¬ 
ficients become mass-dependent: 

Z(gfr,an)^Z(gl,aM,an), c(g§) -> c(g%,aM) . (3.1) 

Here, Z and c are renormalization constants (depending in general on a renormalization scale /i) and 
improvement coefficients, respectively, and we recall that aM is the bare subtracted mass matrix. 
The action in the massive scheme then simply reads 

+ ac sw (gl,aM)xff^(J^ v F^ v \i/ . (3.2) 

This scheme dispenses with the need for the non-perturbative determination of a large number of 
improvement coefficients at the cost of a more complicated renormalization pattern. Note that 
the mass-dependence of the improvement coefficients, Eq. (3.1), is not a necessity for 0(a) im¬ 
provement, as the difference is of 0(a 2 ) in the action. However, we will use this form both for 
consistency and in the hope of reducing the overall 0(a 2 ) cutoff effects. 
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4. Determination of c sw in Nf =3 + 1 

There are two qualifications to be added to the determination of c sw in the massive scheme. 
First, the major effect of the non-vanishing masses on c sw stems from the charm quark. In order to 
substantially reduce the effort, it is therefore reasonable to treat the up, down, and strange quark as 
mass-degenerate light quarks with mass am q j, i.e. to determine c sw in Nf = 3 + 1. 

Second, it is not practicable to vary all the hare parameters (gQ.am q j.am qx ) and find an in¬ 
terpolating formula for c sw in this 3-dimensional space. Hence, we will fix the masses to approxi¬ 
mately physical values, indicated by an asterisk. These two approximations, 

,2 N/ = 3+l , 2 \ fixren. / 9 + * \ . . ,. 

c sv/ (g 0 ,aM) -> c sw (g 0 ,am q j,am qjC ) - > c sw (g 0 ,am ,,am ) , (4.1) 

masses 

together explicitly introduce 0(a 2 ) effects proportional to 

m q . u ~ m q t , m qJ - m q h m q . s - m* ql . m qx - m q c , (4.2) 

which modify the already existing cutoff effects of this order. These modifications are small if the 
mass differences are small. The relevant point obviously is to ensure that the last difference is not 
too large. 

4.1 Line of constant physics 

The improvement condition on c sw will be imposed with the help of Schrodinger functional 
correlation functions (sec. 4.2). The relevant physical scales of the Nf = 3+1 theory in such a finite 
volume are the size L of the box and the renormalization group invariant masses M\ and M c (not to 
be confused with the mass matrix M ) of the light and charm quark, respectively. The bare masses 
are fixed to m ql and m q c (cf. Eq. (4.1)) by requiring the RGI masses to assume their physical 
values, M* and M*. In order to fix these in physical units, a scale is needed. The obvious choice is 
L, which we keep constant as well, at a value L* to be specified: 

L = L\ LM,= L*Mf, LM C = VM* . (4.3) 

Not only is the approach to keep all physical scales fixed convenient, it also has the advantage of 
avoiding possibly large 0(a) ambiguities in c sw , by forcing these to vanish proportional to a in the 
limit —> 0. A discussion can for instance be found in section 1.2.4.1 of [16]. The conditions (4.3) 
define a line of constant physics (LCP), i.e. they determine the bare parameters (gQ,am ql ,am qc ) 
for a given lattice resolution L/a. 

In order to realize the LCP on the lattice, it is reformulated in terms of three more easily 
accessible quantities, 


4>i = g 2 GF (L) 

4>2 = LT u d 

< f > 3 = l - ^r MC - -. 
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(4.5) 

(4.6) 
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Here, is the gradient flow coupling [17] as defined in [18] for Schrodinger functional boundary 
conditions with c = 0.3 by means of the Wilson flow. The functions 

r, 7 = -5 0 log (f'hix o)) \ X0=T/2 , (4.7) 

are effective pseudoscalar meson masses with quark flavors i,j defined in terms of the Schrodinger 
functional correlation function /a (see e.g. Eq. (6.11) in [3] for its definition). All <*>, are associated 
with T = L and vanishing boundary gauge fields. While each of them is a function of all three 
physical scales (or equivalently all bare parameters), they were chosen such that each of them 
dominantly depends on a single parameter. In terms of the <*>,, Eq. (4.3) gets replaced by 

<*>; = <*>*, i = 1,2,3 , (4.8) 

where the (continuum) values <*>* correspond to the scales L*. Mf and Mf. In order to do the 
translation between the two sets of quantities, one needs scale setting, mass renormalization and 
the input parameters Mf and Mf. Since this is not available in Nf = 3 + 1, we resort to Nf = 2 
results. While this alters the particular LCP that is established, it again only concerns the explicitly 
introduced 0(a 2 ) effects in (4.2) via the values with asterisk. In Nf = 2, the setting of the scale is 
available from [19], and the hopping parameters which correspond to the input parameters Mf [19] 
and Mf [20] can be derived under use of mass renormalization constants [19, 21], the improvement 
coefficient b m [21] and the critical hopping parameter fC cnL [22]. We then reused the ensembles 
produced in the course of [22] with 

L* « 0.8 fm (4.9) 

and L*/a = 12,16,20,24,32 in order to evaluate <*>/(a/L*). The results which we employ to define 
the LCP in its form (4.8) are the continuum extrapolated values. We find 

<I>t = 7.31, <!>} = 0.59, <*>5 = 5.96, (4.10) 

fromL* ps 0.8 fm, Mf «M, (iV/=2) /3 = 138/3 MeV [19], and Mf M^ Nf=2) = 1.51 GeV [20]. 

4.2 Improvement condition 

Along the LCP, we impose the improvement condition for c sw in the Schrodinger functional. 
As this follows a standard procedure, we will not cover it here and refer to the original literature [3] 
instead. Concerning the choice of Schrodinger functional parameters, we follow [10], with only a 
few exceptions. First, with respect to the LCP, T = L* is kept, while L is reduced by a factor 2, 
L = Lf / 2. This in particular implies that we use an even number of lattice sites in the temporal 
direction, in contrast to [10]. Furthermore, the improvement condition in our case is 

AM ud = 0. (4.11) 

The quantity A M ud corresponds to AM (cf. Eq. (3.10) in [10]) in the previous determinations of 
c sw at vanishing quark masses. However, in contrast to these, we have to specify the quark flavors 
(u,d) which enter, due to the different masses involved in our case. We thus use the PCAC relation 
in the light sector, as we are mainly interested in improving light quark physics. Moreover, we do 
not use a tree-level value AM* 0 ' on the right hand side of (4.11), since this is not needed when one 
performs the improvement on a line of constant physics, L = L k . 
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C SW 

So 

Kl 

K c 

4>i 

<j>2 

^3 

L*AM ud 

1.6 

1.7188 

0.1403 

0.1222 

7.10(4) 

0.56(2) 

5.94(1) 

0.133(7) 

1.9 

1.7847 

0.1374 

0.1203 

7.18(2) 

0.59(2) 

5.96(1) 

0.046(8) 

1.95 

1.7969 

0.1369 

0.1199 

7.21(4) 

0.60(2) 

6.00(1) 

0.040(8) 

2.1 

1.8340 

0.1356 

0.1190 

7.36(5) 

0.61(3) 

6.02(1) 

0.013(8) 

2.2 

1.8534 

0.1345 

0.1182 

7.17(5) 

0.60(3) 

6.03(1) 

0.002(8) 

2.3 

1.8825 

0.13369 

0.1180 

7.33(6) 

0.60(3) 

5.96(1) 

0.006(10) 

2.4 

1.9069 

0.13278 

0.1175 

7.34(6) 

0.58(4) 

5.94(1) 

-0.018(9) 


Table 1: Preliminary results for L/a = 8. At the bare parameters given in columns 2-4, the values <f>, 
sufficiently close to the LCP were found in T = L = L* simulations with vanishing boundary fields. At the 
very same bare parameters, the quantity in the last column was obtained in simulations with half the spatial 
size and non-vanishing boundary fields. The simulations were performed with the HMC of the openQCD 
package [23] and error estimates have been obtained under use of the T-method [24]. 




Figure 1: L*AM ud (left panel) and (right panel) for L/a = 8 as a function of c sw and their linear interpo¬ 
lation (green) to the improvement point (red). The data point at c sw = 1.6 is excluded from the fits as it is 
clearly outside the range of linear behaviour. 


5. Results for L/a = 8 

Preliminary results for L/a = 8 are listed in Tab. 1. The bare parameters were tuned such that 
the <t>, satisfy the LCP conditions (4.8) with allowed deviations of up to 5% for < Pi, < t >3 and 10% for 
<t> 2 . The results for L*AM ud as a function of c sw are interpolated to the point where the improvement 
condition (4.11) is fulfilled, as shown in the left panel of Fig. 1. A subsequent interpolation of g q 
to the previously obtained c sw gives the corresponding g^, see the right panel of Fig. 1. The results 
of this procedure are 

c sw = 2.26(4), gg = 1.87(1). (5.1) 


6. Summary and outlook 

We have developed a strategy to determine c sw in Nf = 3 + 1 with a massive charm quark. In 
the framework of a massive renormalization and improvement scheme, a line of constant physics 
is set up which aims at keeping in particular the charm quark mass fixed to an approximately 
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physical value in order to avoid large mass-dependent cutoff effects. The first point for c sw as a 
function of gfi on the coarsest lattice with L/a = 8, a « O.lfm, has been computed. Simulations 
on finer lattices, which correspond to smaller g^ and will enable contact with perturbation theory, 
arc on the way. So far, the main effort was in tuning to the proper line of constant physics, but this 
step will become easier due to the gained knowledge. The invested effort is also valuable as the 
approximate knowledge of K c (g q) will help to plan subsequent simulations in the improved theory 
with a massive charm. 
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